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Abstract. Structure and dynamics of complex networks usually deal with degree 
' (— i ', distributions, clustering, shortest path lengths and other graph properties. Although 

O ' these concepts have been analysed for graphs on abstract spaces, many networks 

a^ ■ happen to be embedded in a metric arrangement, where the geographic distance 

! between vertices plays a crucial role. The present work proposes a model for 

growing network that takes into account the geographic distance between vertices: the 
probability that they are connected is higher if they are located nearer than farther. 
In this framework, the mean degree of vertices, degree distribution and shortest path 
Cd ' length between two randomly chosen vertices are analysed. 
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1. Introduction 

For some decades the paradigm of network was represented by the random graphs[lj, 
which modelled graphs with a huge number of vertices and edges as a stochastic process. 
In the well - known Erdos - Renyi model with A^ vertices, the occurrence of an edge, from 
the (2) possibilities, is determined independently with some fixed probability. Despite 
the fact that the introduction of probability concepts in graph theory represented 
a progress in many areas - even besides mathematics -, it failed to describe many 
properties of "real networks" like Internet, air traffic system, social links between people, 
and so on. However, the 1990s witnessed many important progress in attempt to 
examine these "real networks". To cite some of them, the small - world network of 
Watts and Strogatz[2] (see also [3]) is a model that can interpolate, by tuning a single 
parameter, a random graph and a regular lattice. This model succeeded in characterize 
an interval for this parameter where the small - world property is displayed. Barabasi 
and Albert, on the other hand, have proposed the growing network with preferential 
linking[3], which succeeded in display the scale - free properties, present in a myriad 
variety of "real networks". These models have been intensively analysed, with many 
progress. 

Nevertheless, the networks studied by the above authors lied on an abstract space, 
where the geometric structure on which the vertices are placed is not important. This 
picture revealed to be favourable for networks that can be analysed, for instance, by 
tools from spin glass theory or spin systems with finite connectivity that can be studied 
with similar techniques [5|, [6l [71 E]. Many "real networks", however, have an additional 
feature: they are embedded in some geographic structure, which not only allows one 
to distinguish "near" and "far" vertices from some point, but they may determine the 
dynamics of the evolution of networks. For instance, many cities grow starting from a 
"initial" region and expand away from this center gradually; as another example, it is 
common to see people making friendships with those that live nearer than farther. This 
theme - the infiuence of the geography on networks - has attracted attention, and there 
are interesting researches being done in the context of computer science [H [TOl [11], [12], 
with recent progress also in the physical literature [T3| [T^. 

In the literature, networks on Euclidean structures have usually other properties 
attached together p^ [T6| [T7] . For instance, in [16], a random growing network was 
considered on a two-dimensional Euclidean structure together with preferential linking; 
the authors showed, by numerical methods, that the strength of the interaction (that 
depends on the Euclidean distance between them) between vertices can change the 
profile of the degree distribution from a power-law to a stretched exponential (similar 
results were found by |17]), and the small- world property is preserved for any strength 
of this interaction [T7]. A related problem was studied in [18], where the crossover 
between scale-free and spatial networks was considered. The authors of [19] have 
succeeded in constructing a non-growing geographical small-world network with scale- 
free distribution. There are also works that investigated weight properties (differences 



Solvable Metric Growing Networks 3 

in the link "intensities" ) [2U| \W\ on models related to the one studied here. 

Nevertheless, it is important to be aware of the real importance of the metric 
structure on the network. Therefore, a graph which has metric properties only is an 
important model since it concentrates on the influence of geographic structure on the 
network and leads to a better notion of what the real range of a metric network is. For 
instance, scale-free properties is not expected for the model proposed in this work since 
it is purely metric and no preferential attachment is invoked. Searching some general 
features of a network on a metric structure is the main purpose of this work. 

The present work intends to contribute in this direction by introducing a growing 
network that lies on a space with a metric arrangement, where geographic distances 
between vertices are present and do influence on graph properties of the system. The 
model proposed in this work allows one to define vertices that are "far" or "close" 
from a fixed vertex from the standpoint of their geographic arrangement - like what 
happens with a Euclidean space like M^. To stress the metric structure of the model, 
a power - law form, with exponent a, was chosen to describe the probability that two 
vertices, separated by a metric distance r, are connected (the probabihty is proportional 
to r~°). The motivation to choose this particular form is twofold: firstly, the power - 
law form can tune, with a single paramenter (a), the strength of interaction between 
vertices (low values of a makes the network close to a random graph, while large values 
of a makes the network close to a regular lattice structure). In the literature, the 
exponential distribution has already been tested by Waxman|21j (moreover, the random 
network topology generator of Waxman yields a graph that is quite different from the one 
proposed here, where the distance between vertices are "carefully" chosen). Moreover, 
power - law form is interesting since it yielded many non - trivial results to problems 
that are similar to ones that are discussed in this work[9l |10]. The model proposed here 
is manageable to analytic treatments in some interesting cases, and an analysis of the 
mean degree is carried out, as well as other graph properties. 

The next two sections present the model and fix notations. The mean degree of the 
network is analysed in sections IV and V. The aim to discuss this quantity is to show 
that the model proposed in this work exhibits a metric structure that shares the same 
notion of distance present in models embedded in a Euclidean space like, say, M^. For 
sake of completeness and to stress the fact that the model does show metric structure, 
the sections VI and VII discuss, respectively, the degree distribution and an estimation 
of the shortest path length between two randomly chosen vertices in the graph. Some 
final comments are made in the section VIII. 

2. Notations 

Let X be a continuous variable that tends to infinity or some limiting value and let g be 
a positive function of x and consider / as any other function (of x). Throughout this 
work, some symbols [22] will be extensively used, which are: 

I) O: If / = 0{g), then |/| < Ag, where A does not depend on x, for all values of x in 
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question. 

II) o: li f = o{g), then f/g^O. 

III) ^ (or ^): y <^ z will always be understood in the sense that y/ z = o(l) if no 
special conditions are specified. 

IV) 0{x,y): Indicates Taax{0{x), 0{y)} {x and y are assumed to be non - negative). 

In this work, a network will be defined on a metric structure, and the distinction 
between "(metric) distance" and "path length" should be clarified. The "(metric) 
distance between two vertices, xi and X2, is the separation between them on the metric 
structure, and no reference is made to the links between vertices. On the other hand, 
the "path length" between xi and X2 is the minimum number of links that joins them. 
In the example given below (Figure 1), the (metric) distance between the vertices Xa 
and Xh is 5, but the path length connecting them is 2. 




Figure 1. Example: the metric distance between eacti adjacent intersection of dashed 
lines is taken to be 1. The bold curves are links between vertices. 



3. Model 

Let be an initial vertex (the "central node" ) from which the network grows toward m 
branches. At each time step, taken as unity, m new vertices are added to the network, 
one for each branch. Each one of the new vertex (there are m of them: one for each 
branch) is added at distance 1 from the previous vertex in the same branch. Moreover, 
the distance between two vertices located at adjacent branches will be taken as being 1 
if they were born at the same time (this situation is valid for all m, except when m = 1). 
Choose now one branch and denote it by 1; the branches 2, 3, ■ ■ ■, m are counted in 
counter clockwise orientation, as shown in Figure 2. 
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Figure 2. Growing network witli m branches at time t. Tlie probability tliat tlie 
vertex born at time s (located at branch 1) has fc links is Pi(fc, s, t). Note that in this 
example, t — "i and s = t — 1 — 2. By the symmetry of the problem, the probability 
P2{k,s,t) that a vertex located at the branch 2 (this vertex is indicated by an arrow 
in the picture) has k links should be the same: pi(k,s,t) ~ p2(k,s,t). Following the 
same argument, one has pi{k, s, t) ~ p2{k,s,t) = ■ ■ ■ = Pm{k, s, t). 



In this work, the probabihty Pi{k,s,t) that a vertex, which was born at time s 
(< t) at branch i, has k hnks at time t plays the major role. However, the symmetry 
of the model shows that the m functions pi{k, s, t), P2{k, s,t), ■ ■ ■ Pm{k, s, t) are equal. 
Thus, for simplicity of notation, this work will deal with just one of them, say pi{k, s, t). 
Now denote pi{k,s,t) by p{k,s,t). The function p{k,s,t) corresponds to a particular 
choice of one vertex (born at time s) among m without losing generality. This selected 
vertex will be considered to be located, by notation, at the "branch 1". This choice 
does not affect the analysis, since the only important quantity involving branches are 
their differences. 

After t (> 1) time steps, a vertex born at time s, with < s < t, will have four 
nearest neighbors which are at distance 1 (except when m = 1 or m = 2): the vertices 
s — 1 and s + 1 (both of them belonging to the same branch), and the vertices that 
were born at the same time at the (two) adjacent branches. This sets the positions that 
vertices will occupy in the network, but it does not say how they link to each other, 
which will be described below. 

Summarizing the model presented above, one has the following notations: 

a) p{k, s, t) denotes the probability that a vertex born at time s has k links at time 
t. This fixes a specific vertex (denoted as s) among m vertices that have also born at 
the same time. 

b) This fixed vertex s will be considered to belong at branch 1, from which the 
branches are counted up to m in counter clockwise orientation. 

c) By "vertex s" one should understand the vertex which was born at time s or, 
equivalently, a vertex which is at distance s from the central node. 
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d) The parameter t refers to the "time" of the network growing process and also 
indicates the distance of the "border" of the network to the central node. 

Now consider a growing network with m branches and fix a vertex s (at branch 1, 
by convention). At time t(> s), each one of the m new vertices that joins the network 
is made to connect to one of the m(t — 1) + 1 old vertices. The probability that a new 
vertex (born at time t), located at branch b, links to the selected vertex s (at branch 1) 
is denoted by wl {s, t). This probability is made to decay algebrically with the distance 
between vertices, and its explicit form is shown during the next sections. The distance is 
measured according to the previously defined network structure: a vertex s (0 < s < t) 
is at distance 1 from its nearest neighbors along its branch and also at distance 1 with 
its nearest neighbors, born at the same time s, located at adjacent branches; see Figure 
2. These latter two neighbors are absent when m = 1 and m = 2 (when there are only 
two vetices that join the network at each time). The power - law form of wf^ {s,t) 
enables one to observe the effects of the range of the interaction between vertices, as 
can be seen below. 

In this work, rewiring and decimation of edges are not allowed. 

4. One branch 

Consider a network which grows along one direction only (as a "queue"), which is the 
case m = 1. Each new vertex, t + 1 (henceforth, by convenience, new vertices will 
born at time t + 1 instead of t), is linked to just one old vertex, say s, with probability 
w[ {s,t + 1), which depends on the distance (along the line the vertices are placed) 
between them. For simplicity, w[ {s,t + 1) will be denoted as w{s,t + 1) in this one 
branch case; the notation wi™ will be useful later, in the m branches case. 



Figure 3. Rule for the link of the {t + 1) born vertex (here, i + 1 = 3). The vertex 3 
links vertex 0, 1 and 2 with probability w(0,3), 1(7(1,3) and w{2,3), respectively. 

Denoting by p{k, s, t) the probability that a vertex s has k links at time t, the 
master equation of the growing process of this network is casted as 
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p{k, s,t + l) = w{s, t + l)p{k — 1, s, t) + w{s, t + l)p{k, s, t) 



(1) 



where w{s,t) := 1 — w{s,t), and is subjected to the condition p{k,s = t,t > 0) = 6k^i, 
which means that a new vertex has, initially, one ( "own" ) link. 

At time t + 1, the probability w{s,t + 1) of a vertex s be connected by the new 
incoming vertex t + 1 is 

r (t+i-s)-" 



w(sA + l] 



EU r-" 







s < t + 1 



s>t + l 



(2) 



In particular, w{s = t,t) = excludes the formation of "tadpoles" (a vertex that links 
itself). In the one branch case, both the s vertex and the new born t + 1 vertex belong 
always to the same branch. This may not be the case when more branches are present. 
Defining 



{k){s,t):=J2kp{k,s,t) 



(3) 



fc=i 



as the mean degree of the vertex s at time t, its equation of moviment can be casted as 

{k){s,t+l)-{k){s,t) = w{s,t + l) (4) 

by using the master equation ([T]). 

Taking into account the condition {k){s = t,t) = 1 (which is derived from 
p{k, s = t,t > 0) = 6k,i), the solution of this equation can be determined as 

( t 

{k){s,t)=l r=s+l . (5) 



The behaviour of {k){s, t) depends on the order of magnitude of both s and t. The 
analysis will focus on the regime 1 -C s, which may, at first, seem to ignore the region of 
the network "close" to the central node; however, for sufficiently large t (in other words, 
for sufficiently large network), the relative position s, initially chosen to be farther from 
the point 0, is translated to a region near to the central node. Therefore, for the purposes 
of this work, selecting s to be much larger than 1 is not a serious restriction. As a remark, 
the regime 1 ^ s ^ t enables one to cast the continuous version of the equation (jl]), 
although this procedure of transforming a difference equation to a differential equation 
is not necessary, since the solution ([5]) could be obtained. 

Before presenting the results, equipping this m = 1 model with preferential linking 
term (which makes highly connected vertices more likely to receive more links), leads 
to an already studied work[ 
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In the regime 1 ^ s ^ t, which leads the relative position s on the network closer 
to the central node for sufficiently large t, one has 

'^ (l-a)ln(-J +0(1) , 0<a<l 



{k){s,t) = { 



ln('|^)[l + o(l)] , a = l ■ (6) 



The mean degree in a more general regime, 1 ^ s < t, is presented in the Appendix 
lA.ll For a < 1, the mean degree is always increasing with the size of the network. 
Nevertheless, when a becomes larger than 1, the contribution to the degree of a vertex 
comes mainly due to the next vertex only. The calculations have assumed that s <^ t in 
the sense that ln(t/s) 3> 1 and ln(lnt/lns) ^ 1 for < a < 1 and a = 1, respectively 
(theses points are detailed in the Appendix A.l). 

5. m branches 

The general case, where the network grows toward m {> 1) directions, shows that the 
number of branches plays a major role in the behaviour of the mean degree of a vertex. 
The probability that a new born vertex, at branch b (and time t + 1), links to an 
old one, born at time s (and located, by convention, at branch 1), is given by 

''^+ ^ N{t + l,m,a) ' ^^^ 



w 



(m), 



where 



i— 1 m 



N{t, m, a) := t"° + J] Zl ^^"^ ^"' ^' ^'>~" (^) 

u=l h=l 

is the normalization (see Appendix IA.2I) and 

t/^™) (s, t,b) ■.= t-s + min {2s, min{6 - 1, m - 6 + 1}} 

(9) 

is the distance between these two vertices. This is a natural generalization of the previous 
probability ([2]). Here, the distance should take into account the fact that the minimum 
path (over the metric structure) between a new born vertex and some old vertex may 
cross the central node (this possibility is not seem for m < 5). This explains the term 
min {2s, min{5 — l,m — b+ 1}}, where 2s stands for a path (over the metric structure) 
that passes along the central node, and min{6 — 1,?ti — 6 + 1} is appliable when the 
shortest "metric" path realizes contourning the central node. For the model with m 
branches, the master equation is written as 

m m 

p(k, s,t + l)= p(k - m, s, t) Yl w\""^ (s, t + 1) + p(k, s, t) JJ wf"'' (s, t + 1) + 
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m—l 



+ J2p{k-b,s,t) Yl w^^\s,t+l) 



■w 



Ml 

Qb ' 



sA + l) 



n 



w 



(m) 



S,t + 1) 



(10) 



6=1 



qi<-<qb 






where, in the right side, the first term represents the case when all the (m) new vertices 
links to a specific one s and the second term is when none of the m new vertices links 
s. The third term is the sum of the remaining intermediate cases. 

As in the previous section, one is interested in the mean degree of a vertex s at time 
t, namely {k){s,t). From the definition ([3]) and the master equation flTO|) . it is possible 
to show, after a lengthy arithmetic manipulation, that (see Appendix IA.3J) 



{k){s,t + l) = {k){s,t) + Y,w'r\s,t+l). 
The solution of the difference equation (TTTl) is 

/" m t 

1+E E^^^(^'^) 



:iii 



{k){s,t) = { 



6=1 r=s+l 



S <t 



t 



(12) 



To stress the influence of the number of branches on the graph property of the 
network, the mean degree will be calculated in the regimes 1 < m ^ s ^ t, where the 
number of branches is small, and 1 -C s <C t -C m, where m is large. Actually, in the 
former case (1 < m <C s <^ t), one has the same results for the m = 1 model: 

^ /A 

a = 



{k){s,t) = { 



\n{-]+0{l) 



;i-a)ln( - ) +0(1) , 0< a< 1 



(13) 



ln(;^) [1 + 0(1)] 



2 + Ois~\T 



1 ^l-a\ 



a = 1 



a > 1 
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10 



On the other hand, if 1 ^ s ^ t ^ m (in the sense also that In (t/s) ^ 1), one has 



' lnfi]+Ois 



{k){s,t) 



-1^ 



)l-a 



^.„ il+od) 



ln2 Vs 



0{1] 



a = 



< a < 1 



a 



(14) 



^^ln(i)+0(l) , l<a<2 



a 



1 



'^J-(-.}+oi^- 



a>2 and m > t""^ 



a > 2 and m <^ s 



a-l 



2 + 0(msi-°) 

The results for the mean degree in other regimes are presented in the Appendix 
IA.4I for completeness. Note that the asymptotic behaviour of the mean degree for 
large number of branches is distinct from the one branch (or relatively small number of 
branches - in the sense that m -C s) case, showing that m does influence on the graph 
property of the network. Basically, one has {k){s,t) ~ 21-^-1 ^^ (s) ^°^ ^^^ ^ ^°^ large 
m, except when the number of branches is large but not larger than s"~^ (see the a > 2 
case). Actually, this last case shows that despite the fact of m being large, it belongs to 
the case where the interaction parameter a is strong enough to let each vertex to have 
two links only: one from its "own" and the other from the vertex born immediately 
after it in the same branch; in other words, the interaction a overcomes the size of the 
network and each vertex can see its neighborhood only (therefore, when a is sufficiently 
large, the network effectively behaves as a small m case). 

As a last remark, the number of branches became important for non-zero a only, 
which leads the model to "notice" the metric structure of the space. In the particular 
case where s -C t, in the sense that ln(t/s) ^ 1, one has {k){s,t) = ln(t/s) + 0{1). 
These mean degrees recover results from recursive random graphs, where the metric 
structure that the system lies on can be ignored - this is cleary seen from the fact that 



w 



(m) 



s^t + 1) = {mt + 1)"^ for a = 0. 



6. Degree distribution 

Consider first the case m = 1. From the master equation ([T]), it is possible to show that 
(see Appendix IA.5I for technical details) 

fc-i 



p{k,s,t) 



{k-l)\ 



^ w;f^(s,r) 



r=s+l 



exp 



E-l 



(1), 



S,U) 



u=s+l 



[1 + 0(1)] 
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if 

iv'i\s,u) -^1, for M G (s,t] C Z+. (15) 

This condition can be satisfied for < a < 1 if s ^ 1, but it does not hold for larger 
values of a, except if a ~ 1"*". The form of the equation for p{k, s, t) given above shows 
a Poisson distribution with mean ^^^g_,_]^ w\ (s, r) = {k){t) — 1 for s < t (see (|5])). This 
means that the probability of a vertex s, at time t, having k links is concentrated on 
the mean {k){t) — 1, and the structure of the network is homogeneous. 

Using the asymptotic form (fT5|) . one can show that the degree distribution P(k) : = 
\imt-^oo P{k,t), with P{k,t) := J2s=iPi^^^^^)/^^ behaves as 

fc-i 



P{k) 



1 n-a\ 

< a < 1 



a \2 — a 

(16) 



(fc-1)! 
For a = 1, it is possible to evaluate the degree distribution at time t as 

P{k, t) = 2-'-^ fcFfe(2, ■ ■ ■ , 2; 3, ■ ■ ■ , 3; In t) 
Int^ 1 (Int)" 

where kPk{-, ■', •) is the hypergeometric function. 

A simple calculation shows that for s ^ 1, one has dp/ da < for sufficiently large 
k. The degree distribution decays more rapidly as a increases. For instance, in the case 
< a < 1, the decayment is exponential, and for a ~ 1+ the degree distribution decays 
more rapidly, as an inverse of a factorial. Moreover, in the limit of a ^ oo, P{k) tends 

to 6k,2- 

For the other case, when the number of branches is large in the sense that m ^ t, 
a heuristic approach will be adopted to estimate the form of the degree distribution. If 
one assumes that the probability that a vertex s has exactly {k){s, t) links at time t - the 
"(5-Ansatz" -, it is possible to show (using results for the mean degree (A;)(s, t)[24l [25] ) 
that P{k) has an exponential distribution 

P{k) ~ ^e-'^/^^-) , (18) 

K[a) 

whenever {k){s, t) ~ K{a) ln(t/s), where K{a) is a factor that depends on a. Assuming 
the "5-Ansatz" , which turns to be reasonable in the regime 1 ^ s ^ t, the exponential 
form of the degree distribution suits the cases where m ^ t ^ s ^ 1 (for any non- 
negative a), except when a is sufficiently large (in the sense that ms^~°' <^ 1), when 
P{k) tends to 6k,2- 

When each vertex has just two links on average due to the high value of a, the 
model indicates that each vertex has, apart from its "own link", another that comes 
from its successor in the same branch. Thus, the network achieves a regular structure. 
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like a regular lattice, where many of the vertices have connection to its nearest neighbors 
only. 

The metric nature of the graph is well characterized in the regime where the number 
of branches is not large like m ^ t. Despite the fact that the probability that decays 
algebrically with distance can tune the network between a random graph type one (low 
values of a, typically < a < 1) and a regular lattice type (large values of a, when 
{k){s,t) ~ 2), if the number of branches is too large, the normalization ([8]) diverges 
with m and w^™ (s,t) becomes "nearly" uniform. This is the origin of the observed 
randomness in the network connections, which reflects in the exponential form of the 
degree distribution ( ITSll . 



7. Shortest path length 

This section will provide an estimation for the mean shortest path length i{t) between 
two randomly chosen vertices, say, Xa and x^- Since they are randomly chosen, their 
distance from the central node will be taken as 0{t). 

In this work, the order of magnitude of the shortest path length between Xa and x^ 
will be considered to be the graph distance (path length) between Xa and the central 
node plus the graph distance between Xh and 0. Denoting the (shortest) path between 
two vertices, Xi and X2 by Xi ^^ X2, the estimation of the shortest path between Xa and 
Xb proposed in this work stems in replacing Xa ■«-»■ x^ by Xq <-*■ <-*■ x^. This estimation is 
expected to be reasonable at least in the case where the number of branches is sufficiently 
large, because then the probability of the paths Xa ^^ and Xb ^^ having intersection 
other than the central node is low (meaning that the shortest path Xa ^^ Xb includes, 
necesarilly, the vertex 0). It will be shown later that if m is sufficiently larger than t, the 
shortest path length is estimated as 0{lnt), which resembles the result from random 
graphs. Note that by construction, there are no two different shortest paths connecting 
two points in the model presented in this work, since a new born vertex links to only 
one old vertex. 

Next, define the function 

r— 1 m r— 1 m 

s=l &=1 s=l 6=1 

which is the measure of the mean (metric) distance covered by a link (starting from a 
point at distance r from the origin) toward the central node. In other words, this is the 
metric length of the projection of a link on some fixed branch (see Figure H]). 

The evaluation of A^ (r) will focus on two distinct limit situations for the number 
of branches, and it will play a major role in the estimation of the shortest path length. 
As mentioned earlier, one should estimate the (order of the) number of paths that 
a randomly chosen vertex Xa takes to achieve the central node (and the same for Xb)] 
furthermore, as stated in the beginning of this section, it will be assumed that Xa = 0{t) 
(the distance between Xa and 0) and by the same argument, one has also Xb = 0(t). 
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Figure 4. Let the path xi ^^ a;2 ^^ 2^3 ^^ ■ • • be the shortest path between xi and 
the central node. The projection of X2 and X3 (and also xi) onto the branch that joins 
xi and the central node are represented by a filled ball in the figure. The quantities 



were not drawn 



A\,"''{xi) and A^ {X2) are also indicated. For sake of clearness, most of the branches 



This means that the path length Xa ^->- ^^ Xb (which is an estimation of the shortest 
path length Xa ^^ Xb) is of the order of the path length Xq <-> in this work. 

Let Xa = Xi <-H> X2 <-^ ■ ■ ■ <-* x„ = be the shortest path between Xa = Xi and the 
central node = x^. Using ( TT9l) . one has 

n-l 



E^ 



Ml 



Xi 



distance between Xa and 0, which is 0{t) . 



(20) 



Remembering that n is of the order of the shortest path length i{t) between Xa and Xb 
at time t, the relation (|20|) can be used to estimate i{t). 

The details of the estimation is presented in Appendix I A. 61 For m <^ s, when the 
number of branches is not large, one has 



a 


m 


<a < 1 


C(lnt) 


1 <a< 2 


o(r-i) 


a = 2 


0(t/lnt) 


a>2 


0{t) 



On the other hand, if m ^ t, one has 



a 


m, t 


£(t) 


0<a <2 


— 


O(lnt) 


a > 2, 


mt^-" > 1 


O(lnt) 


2 < a < 3, 


mt^-" < 1 


O(r-Vm) 


a = 3, 


mt^-" < 1 


C(tV(m In t)) 


a > 3, 


mt^-" < 1 < mt2-° 


C(r-Vm) 


a > 3, 


mt^-" < 1 


0(t) 
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The above results indicate that whenever the number of branches is too large (in 
the sense of mt^~°' ^1), the model behaves as a random graph, leading to a shortest 
path length between two randomly chosen vertex that increases logarithmically with t. 

8. Conclusions 

The present work has introduced a network growing toward m branches embedded on a 
metric structure and analysed some graph properties. The mean degree of a fixed vertex 
decreases as a increases, indicating, as expected, that the strength of the "interaction" 
confines the vertices to link to other vertices that are located nearer than farther. To 
ensure the fact that the number m of branches does influence on the graph properties like 
mean degree, two opposite conditions were then considered. For small m, the results 
were similar to the one branch case, while for large m, the mean degree was larger. 
In any case, for a sufficiently large, one has {k){s,t) -^ 2, no matter the number of 
branches. The degree distribution was evaluated for < a < 1 and a ~ 1"^, and it 
shows, as expected, that they decay as rapidly as the magnitude of a increases. Finally, 
the shortest path length between two randomly chosen vertices was estimated as a 
function of t. The metric structure of the model is better displayed if the number of 
branches are not overwhelmingy large; otherwise, the model behaves as a random graph, 
and this fact is supported by the results obtained for the degree distribution, shortest 
path length and mean degree; this last one grows as ln(t/s), when s <C t, for any a for 
sufficiently large m - see (1141) -, as it is expected |25j. The model presented in this work 
was successful in displaying properties that are common to networks that are embedded 
in a Euclidean space (like M^) and is manageable to analytic treatment. 
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A. Appendix 

This appendix is devoted to present some technical observations and list results omitted 
in the main text for clearness. 

The present work has used extensively the fact that 

-^^ + C(a) + C(r") , 0<a<l 
1 — a 

lnt + 7 + C(ri) , a = l , (21) 



t-i 



C(a) - ^^ + 0(t-") , a>l 
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where ( is the zeta function of Riemann and 7 (= 0.5772...) is the Euler - Mascheroni 
constant. As a remark, ({a) > for a > 1 and C{a) < for a < 1. 



A.l. One branch case 
For 1 ^ s < t, one has 



l + (l-a) 



dr 



r (1 — r) 



s/t 



+ C>(s"-^) 



< a < 1 



{k){s,t) = { 



1 + 



dr 



(r — s) Inr 



+ C(ln-^s) 



(22) 



s+l 



^-^,p-'-Ois^-) 



a = 1 



a > 1 



where t := t — s. 

The case a = 1 may deserve a short comment. Consider the case a = 1, m = 1 
and 1 <^ s -C t such that t — s :^ s'' for some positive 6 = 0{1). Using (12T|) . the mean 
degree ([5]) is 



{k){s,t) = i+ J2 



r — s] 



r=s+l 



lnr + 0(l) 



t-s 



dy 



t-s 



y\n{y + s] 



+ o\i, 



dy 



yln {y + s) 



(23) 



However, since 



t-s 



dy 



yln (y + s) 



< 



dz 



t-s 



+ 



dz 



zln s J zln z 



In(t-s) 

0(ln-'s)+ j ^ = 0(1) 



(24) 



one is left with 



(fc)(.,t) 



dy 
ylniy + s] 



+ 0{l) 
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t-s 



dy 



+ 



y\n{y + s) J y In {y + s] 

1 s'5+l 



dy 



+ 0{1) 



where the first term (last hne) on the right hand is bounded {6 = 0{1)): 



„*+i 



dy 



y\n{y + s) 



j+i 



Ins y y 

1 



(25) 



(26) 



Thus, 



t-s 



{k){s,t) 



dy 



t/s 



y\n{y + s) 
dz 



+ 0{l] 



t/s 



1+s 



z\u.[sz) 

dz 
z\n{sz) - 



1 



z-l 



0{l) 



l + 0{s-') +0(1) 



(27) 



and the desired result follows. 

The estimation for a = 1 assumed that t ^ s^ for some positive 5 = 0{1). 

A. 2. On the normalization ^ 

The (asymptotic form of) normalization factor is evaluated by considering the cases m 
odd and m even separately. It is shown that both have the same general formula. 
For m odd, the normalization is 



r-l 



Tn — 1 

r—1 2 



Noddir,m,a) = r " + ^ (r - s) " + 2^ ^ r - s + min{2s, j} 

S = l S = l J = l 



f28) 



where the first term is the distance between a r born vertex to the central node; the 
second term counts the distance between this new vertex to the vertices located in the 
same branch, and the last term the remaining ones. 
If m is even, one has 

r— 1 7 — 1 



Neven{r,m,a) = r °' + y^{r — s) " + / ^ | r — g + min{2s, m/2} 

s=l s=l 



+ 



+2 y^ y^ \r — s + min{2s, j } 

s=l j=l 



(29) 
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where the meaning of each term is similar to the odd case. The third term is the 
contribution due to the branch that diametrically opposes branch 1, which is absent if 
m is odd. 



For 1 < m '^r, one has 



mr 



l-Q 



1 — a 



+ C(m,a) + 0(r^°) , 0<a<l 



Nr odd |(r, m,a) = < 

\ even J 



where C{m,a) is defined in flTTI) . 
For m ^ r, one has 

A^r odd \{r,m, a) = 



mlnr + C{m,l) + 0{r ^) 



a = 1 



, (30) 



C{m, a) 



mr 



l~a 



a — 1 



+ 0(r-") , a > 1 



If odd T 
\ e^en J 



)1 — a 



1 — a 



mr 



i-« /■l + 2-n__„ , 8[(l + a)2-"-l] . „ 



mr + 



[1 -a) (2 -a) 



-r'""+ 



— o— 1 ^l—a\ 



+0{mr " ,r 



mln2 — —mr ^ + A (2\n2 — I) r + 0(\nr,mr ^) 



< a < 1 



a = 1 



1-2 



l-a 



a — 1 
+0{mr-'^'\ 1) 



i^a /l + 2"°\ _„ 8[l-(l + a)2-] 2_„, 



(a-l)(2-a) 



(31) 
1 < a <2 



1 5 

-mr'^ mr"^ + 2 Inr + 4 + 27- C(2) -6 In 2 + C>(mr~^r~^) , a = 2 



1 _ 2^^"\ /I -I- 2~"" 

mr^-'' - I —- 1 mr-" + 2C{a - 1) - ((«) + 



a — 1 



8[(l + a)2 - l]^2^,_^^^^^_„_i^^i_„^ 



a >2 



(a - 1) (a - 2) 

74.5*. On i/ie m branch master equation 

After multiplying both sides of the master equation (TTOl) by A; and summing them, one 

gets 



{k) (s, t + 1) = A(A;)(s, t) + B{s, t + 1) 



(32) 
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Since one can show that 



A := l[w\"'\s,t + 1) + Hwf )(.,t + 1) + 



+ 



m— 1 



Z^ z^ 

b=l qi<---<qi 



W 



(m)| 
91 ' 



S,t+1)- 



■W 



Ml 

Qb 



s,t + l] 



n ^"^^^,^+1) 



i=i 

j¥=gi,-,gb 



(33) 



and 

5(s,t+l) := 



m— 1 



771 



J=l 6=1 qi<---<g6 



'Krn^,^+i)x 



X 



n 



mr'(s,(+l) 



j=l 

j¥=Cll,-,<lb 



E 



w;f^(s,t + l) 



(34) 



for all positive integer m, the results follows. The formula (1331) can be easily obtained if 



one notices that ^ = H 



m 



w- '(s,t+i)+w^r\s,t+i: 



(m) . 



. On the other hand, although 



the deduction of the formula flMl) is straighforward, it is tedious. It may be of some help 
to introduce the quantity r^ (s,t) := Wq {s,t)/w^I^\s,t), which is well defined (since 



Wr. 



{s,t) does not go to zero), and deduce the formula by induction by working out the 



expression 



B{s,t + l) 



m 



n 



w 



(M) 



(s,t + l) + 



n 



wf'\s,t + l] 



X 



(35) 



m—l 



x2_^b 



5^ r(f)(.,t + l)...r(f)(.,t + l). (36) 

6=1 qi<---<qb 

Note that the superscript of the probability Wj {s, t+ 1) was replaced by a fixed number 
M. The induction argument will be applied to the variable m only, and it can be shown 
that it is valid for any M - in particular, in the case of interest M = m. The reason 
for make M independent of m is twofold. First, the expression (IMl) is valid for any 



w 



(m). 



(m) 



(m). 



s,t + t) and w, (s, t+ 1) that satisfies w:- {s,t+l)+w) (s, t + 1) = 1 for each j 



tM/ 



and does not depend on the particular choice of wf" [s, t + 1) and wf^ (s, t + 1) like ([7]) 



(m). 



SM) 



(M), 



to define r\ {s,t + 1), one should impose the additional condition Wj {s,t + 1) 7^ 0. 
Secondly, without making M independent of m (in other words, if one keeps m instead of 
M), the induction argument became extremely tedious to be applied, since the function 



w, 



(m) 



s,t + 1), from the induction hypothesis, is different from the w 



(m+l) 



s,t + l) 



see 



the definitions (171) and (IHl) 
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Before analysing the mean degree for many cases not presented in the main text, it is 
convenient to define 

■m — l 



mCia) -2J2Y. 



u 



m 



C{m, a) := < 



6=1 M=l 

m/2 f-1 b 

c(«)-E-""-2EE-" 

u=l b=l u=l 



m odd 



m even 



> a^l 



.(37) 



m — l 

~^2~ b 



m7-2^^M ^ 



, m odd 



6=1 M=l 



> a = 1 



m even 



m/2 f-1 b 

m7 - ^ M"^ - 2 ^ ^ n"^ 

M = l 6=1 u=l 

This work has concentrated in two hmiting regimes to stress the role played by m. 
The results for "small" number of branches {1 < m '^ s < t) was compared to the other 
limit case, where m is "large" (1 <^ s < t <^ m). In both cases, s is chosen to be much 
larger than 1 and in this appendix, t(> s) is taken free, although simple analytic results 
are provided if s <^ t, as shown in the section [51 

For 1 < m -C s < t, one has 



{k){s,t) 
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1 
1 + (1 - a) / ,/^ ,^ + 0(s"-^) , < a < 1 



1 + 



s/t 



dy 



s+l 



{y- s)\ny 



+ C(ln-i s), a = l 



m({a) 1 

C{m, a) C{m, a) 



m — 1 

~2~ b 



m—l 



25:5:"-"+ E"-" + 2E E 



M 



6=1 M=l 



M=T+1 



b=l u = T + b+l 



(38) 



+ 



+0(s^-"), modd, a>l 



1 + 



m^(a;) 1 

C(m, a) C(m, a) 



6=1 n=l 



n=l 



6=1 u=T+6+l 



+ E"-°+ E 



M 



n=T+l 



u=T+m/2+l 



+ C(s 



l-a\ 



m even , a > 1 



where r := t — s. The calculations can not be always computed directly in the a > 1 
case, and the mean degree is estimated by evaluating bounds. However, in many cases 
the lower and upper bounds have the same order of magnitude with the same leading 
term. The results in ( !38l) show that the behaviour of the mean degree for 1 < m <^ s < t 
is similar to the one branch case. 

On the other hand, for 1 ^ s < t ^ m, one has 



{k){s,t) 
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In ( ^ ) + 0{s~') 



1 + 



1 — a 



1 



n-a 



dy 



U yii + yf 

s/t 



+ 0{s^\t/m) 



a = 



< a < 1 



ln2 V 2s 



+ 0[s ^,t/m,{s/m)lnt,(T/m)\ns 

, 1 < a < 2 



a 



m 



1 



a — 1 



dy 



m s 



l-2i""7 vil + yY 

s/t 



+ < o 



o 



tint 1 



m 



a 



o-l 



a >2 



where it was assumed that slnt -C m and rlnt <^ m for a = 1; tint <^ m for a = 2; 
and t°~^ <^ m for a > 2. 

For 1 ^ s < t ^ 771, the interplay between the free parameters should be examined, 
specially the competition between the strength of "interaction" a and the number 
of branches m. The present work has considered the limiting cases t°^^ ^ m and 
^a-i ^ fYi iq show their differences. Since m ^ t, note that the former condition 
("large m" - in contrast to the later, s°~^ ^ rra, the "large a") is always satisfied for 
< a < 2. Finally, the parameters mt^~°' and ms"'^ dictates the behaviour of the 
mean degree in the case a > 2. If t"~^ <^ m, one has the situation shown in (l39l) . where 
for sufficiently large m (such that m ^ t"^^), the mean degree grows with the size of 
network. On the other hand, if s"~^ ^ m, or, equivalently, the a is sufficiently strong, 
one has fl38l) . where the mean degree behaves as in the case m = 1. 

A. 5. On the degree distribution (m = 1) 
Performing the zeta transform 

oo 

P{K,s,t):=Y,K'p{k,s,t) (40) 

fe=i 

on the master equation ([1]), one can show that 

t 

= kY[[1 + {K-1)w{s,u + 1)] , 



(41) 
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since p{K, s, m) 7^ for any u > s. From the condition p{k, s, s) = Sk^i (or p{K, s, s) 
K), one has 

fc-i 



p(/c,S,t + l) = y^;— — ^1 1 X 



(k- 



X 



q=0 

dK'^-^-i 



k-1 



exp 



{K -l)^w{s,u + l) 



K=0. 



T\a,K) 



where 



' O (g!s2("-i)) , < a < 1 



T\a,K) 



5q,0 + < 



K=0 



O 



In^s 



a = 1 



a ^ V 



K=0 



,(42) 



(43) 



for sufficielntly large s. Therefore, for large s (say, s ^ e^ ), one has 



p{k,s,t) 



-m{s,t)-i] 



k-l 



mis,t)-ir-'[i+oii)] 



where the main contribution comes from the q = term (in the sum 
is given by ([5]). 



(44) 
and {k){s,t) 



A. 6. Shortest path length 

Firstly, the asymptotic form of As (r), which is closely related to the estimation of the 
shortest path length, will be listed. The definition fl371) will be used below. 
For 1 < 771 -C r, one has 

~ " ' r + 0{r'' 



and 



A(™)| 



2-a 



— + 
mr 



In^r 



< a < 1 



a = 1 



r = < 



(45) 



mr 



2-a 



C{m, a) (2 — a) 



+ 0{r-^-"^, 1) , 1 < a < 2 



m 



[ C(777,2) 



lnr + 0(l) 



a = 2 



C{m, a — 1) 
C{m, a) 



+ 
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1 v:^ fm^ - 1 



C(m, a) ^-^ 

^ ' M = l 

C(m,a) \ 4 



1 ^" 
t V 

C{m,a) ^ 



+ u-u'^ U""- 



m odd 



m^ — 2m 
4 






m/2 

M=l 



+ m-m2 )«-"+ 



(^(a) m fm 



> +0{r 



2-a\ 



C{m,a) 2 V2 



m even 



(46) 



for a > 2. Note that As (r) becomes smaller as it approaches the central node. 

On the other hand, for 1 <^ r <^ m, the parameter mr^~°' and mr'^~°' play a major 
rule to the asymptotic behaviour of As (r) . Firstly, one can show that 



AM I 



fa + 2 



2-a 



3)r 



, C 1 

(2 - a) (21-° - 1) V m 



2-7^ V + c(l,— 
ln 2 / \ m 



< a < 1 



a = 1 



r) = < 



-^ '-— + {1,— ] , l<a<2 

2 - a) 21-° - 1 V ^ 



(47) 



(2-21n2)r + C 1, 



r^lnr 



m 



a = 2 



(4^ 



where a stronger condition, m ^ rlnr, was assumed for a = 2. For a > 2 and 
mr^^°' ^ 1, one has 

r (a + 22-°-3)r 



AM| 



r = < 



X > C* 1,— 
(2 - a) (21-° - 1) V m 



2 / r^ 

Trr + O 1,- 

3 \ m 



I (2-a) (21-° - 1) V ' ^ 



2 < a <3 



, a = 3 



a > 3 



(49) 



(50) 

In this regime, one sees that the leading term of A^ (r) (for large r) is proportional to 
r. Note that the condition mri-° ^ 1 is automatically satisfied for < a < 2 (with 
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m ^ r). A different asymptotic behaviour is seen for mr^ "^l;if2<a<3, one gets 
( (a + 22-" -3) mr^-'' 



AM( 



r = < 



(2-a)(l-a)[2C(a-l)-C(a)] 
1 mr~^ 



+ 



2 < a < 3 



(51) 



+ o{ 



I 4 [2C(2) - C(3)] 
and for a > 3, one has 

f (a + 22-° - 3) mr^-'' 



2 — S — 1 

m r ,mr 



, Inr) , a = 3 



AM, 



+ 



r = < 



(2-a)(l-a)[2C(a-l)-C(a)] 

+C (mV^-^^, 1) , mr^-" < 1 < mr^-" 



(52) 



C(«-2) 



+ 0{mr 



2-a\ 



mr " ^ mr " ^ 1 



(m) 



I 2C(a-l)-C(«) 
From the resuhs, fHSj) to (152|) above, one sees that the asymptotic behaviour of 



As (r) is one of the following form: 



CiT 

r 
Inr 



A^HO~<; c^r^ 



(53) 



Cslnr 



C4 



where Ci (0 < Ci < 1), C2, C3, C4 and cj do not depend on r. For instance, in the regime 
1 < m <C r and for 1 < a < 2, one has C2 = m/[C{m, a) (2 — a)] and cj = 2 — a. 

Each asymptotic form for As (r) leads to a different shortest path length i. 

The shortest path length between two randomly chosen vertices, Xa and Xb, will be 

estimated to be of the order of the path length between Xa and 0, as stated in section 

[71 Symbolically, this implies 

i 

y A^J"'{xi) ~ distance from Xa = Xi to the central node, which is 0{t) , 



j=i 



(m). 



where X2 is closer to the central node than xi = Xq by a distance As ' [xi) on average, 
X3 is closer to the origin than X2 by a distance As (^2) on average, and so on. 

Since xi is at distance 0{t) from the origin on average, it means that if As (^) = 
C(l) (or Ai™)(r) ~ C4), then 

i{t) = 0{t) , when A(™)(r) ~ C4 . (54) 
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On the other hand, suppose that As (r) ~ cir, where ci G (0, 1) is a constant. In this 
case, one has the following sequence: 

X2 = Xi- A^"") (Xi) = Xi- CiXi = {1 - Ci)xi 
X3 = X2- C1X2 = (1 - Ci) X2 = (1 - Ci) Xi 
Xi = X3 - C1X3 = {1 - Cif Xi 

Xn = (1 -Ci)"'~^Xi (55) 

From (!55|) . one has 

n = l-\ y ^ , Ci G (0, 1) . (56) 



In 



1 



l-c 



Remembering that Xi ~ C^(^) and if n is of the order of £, then one has £{t) ~ Int. 
However, this kind of argument is suitable for n not so small (since the estimation of 
As (r) is made for large values of r). Then, to estimate the shortest path length when 
As (r) ~ cir (ci G (0, 1)), n will be taken as being of order of i, which means that x„ 
is close to the central node but not the central node itself; in other words, x„ ~ C'(l) 
and n ~ £. For large values of t, one has 

£(t)~lnt, when A^J^^r) ^ civ for ci G (0, 1) . (57) 

For technical reason, the shortest path length for the other three cases will be 
estimated by using (As (r)) := ZIm^s ('")/^ ~ fAs{u)du/t instead of As (r). 
From this approximation, one has 

£{t) ~ Int, when A^'"^^) ~ 7^ , (58) 

mr 

i{t) ~ t^-^ , when A(™)(r) ~ C2r'' , (59) 



and 
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